Spin Diffusion in Double-Exchange Manganites 
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The theoretical study of spin diffusion in double- exchange magnets by means of dynamical mean- 
field theory is presented. We demonstrate that the spin-diffusion coefficient becomes independent 
of the Hund's coupling Jh in the range of parameters JhS > If > T, IV being the bandwidth, 
relevant to colossal magnetoresistive manganites in the metallic part of their phase diagram. Our 
study reveals a close correspondence as well as some counterintuitive differences between the results 
on Bethe and hypercubic lattices. Our results are in accord with neutron scattering data and with 
previous theoretical work for high temperatures. 
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Spin diffusion dominates the low-w, low-k excitation 
spectrum of a magnet in its paramagnetic state and con- 
tains important information about the spin dynamics. In 
the past, spin diffusion was studied intensively in Heisen- 
berg systems and, recently, has been investigated both 
theoretically and experimentally for strongly-correlated 
itinerant magnets^*^ The current growth of interest in 
spintronics requires understanding how local spins relax 
through their interactions mediated by itinerant charge 
carriers rather than through their direct interactions with 
each other A Among such systems are the colossal mag- 
netoresistive (CMR) manganites which are ferromagnetic 
metals in a large part of their phase diagram^ Recent sys- 
tematic neutron scattering experiments on the ferromag- 
netic CMR materials revealed a peak centered at uj = 
associated with the spin diffusion. This peak was also 
seen below the ordering temperature indicating electronic 
inhomogeneity with regions having lower T c 's. 

In this Letter we present a comprehensive, self- 
consistent, microscopic calculation of spin diffusion ap- 
plying dynamical mean-field theory (DMFT) to the 
double-exchange (DE) model. We demonstrate that the 
spin-diffusion coefficient T> s is related to the local single- 
particle Green's function and can be evaluated as a func- 
tion of doping and temperature. Our results agree quan- 
titatively with neutron scattering data on manganites for 
a range of doping concentrations. Thus, our approach 
creates a framework for the self-consistent study of dif- 
fusive spin dynamics in many real materials, including 
magnetically doped semiconductors. 

Following the general hydrodynamic arguments of 
Ref. 0, we write the generalized susceptibility of a para- 
magnet for low energies and long wavelengths as 



X(q,w) ^ x(q) 



(1) 



— iuj + T> s q 2 ' 

which through the fluctuation-dissipation theorem 
gives the neutron-scattering dynamical structure factor 
5(q,cj) ~ 2[n.B(w) + l]Im%(q,cj), where x(q) is the static 
susceptibility at wave- vector q and ub{u) = [e w / T — 
We take % = Ub = 1 throughout this paper. Further, 
the generalized susceptibility can be related to the spin 



current-current correlation function using the dispersion 
relations 7 " and the continuity equation dS a (r,t)/dt = 
-Vd?(r,t): 



x(q,w) 



where T^f(q,u) = -i J dte^6(t)([jf (q, t), #(q, 0)]) 
is the retarded current-current correlation function®, 
jf (q, t) is the ith component of the spin current for the 
a-spin projection, i = 1 . . . d, d is the dimensionality, 
and a is the lattice constant. We use the isotropy of the 
spins above T c and assume the isotropy of real space to 
suppress the indices in II^ 3 (q, u>) = II(q, u)5 a p5ij . Com- 
bining Eqs. and @ in the q, u> — > limit, we write 
the Einstein relation between the spin-diffusion coeffi- 
cient and the spin conductivity a s (which in general is 
distinct from the particle conductivity) as 
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n(q, w )-n( q ,o) 



(2) 



lim 



lm[n(0,w) 



(3) 



where x = x(Q = 0). These expressions are general and 
do not depend on the microscopic model. 

We now consider the DE model with Hamiltonian 



n 



-t 



H.c. 



2J H S > 



(4) 



where t is the nearest-neighbor kinetic energy, Jh is the 
Hund's coupling between the local Mn 3+ 5 = 3/2 spin 
and the electronic spin s = and <r are the 

Pauli matrices. To describe the multitude of phases in 
manganites requires that the orbital, phonon, or Jahn- 
Teller terms be included in the above model. 9 How- 
ever, the magnetic properties of these materials in the 
metallic part of their phase diagram, such as the mag- 
netic excitation spectrum and the ferromagnetic transi- 
tion temperature^iii^ are quantitatively well described 
by the model in Eq. Q. Therefore, such a model 
must be also capable of describing spin diffusion in these 
systems. 13 Regardless of modifications of the DE model 
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needed to describe a particular real system, elucidating 
the dynamic properties of this basic model of strongly- 
correlated itinerant magnets is an important task on its 
own. 

Within the DE model there is no direct interaction 
between the local spins. Thus, the total on-site spin 
Sj ot — Si + si commutes with the exchange part of Eq. 
l|Hl and the spin current can be expressed in terms of elec- 
tronic operators only: 14 j'f(q) = £ k v l c l..j^s c ^-^/2, 

where v k = V*£k and £k = — 2t cos This is 
simply another way of saying that electrons mediate the 
magnetic relaxation processes in an itinerant system. 

The physical situation relevant to manganites corre- 
sponds to strong Hund's coupling JhS ^> W. Since the 
characteristic relaxation time for the electronic spin is 
short and the spin relaxation is essentially local, pertur- 
bative approaches to the spin diffusion 1 ^ are inapplicable. 
Therefore, we employ DMFT, which takes into account 
the local dynamics in strongly-correlated systems and has 
been successfully applied to a number of problems^Siil 
Using DMFT also simplifies our problem significantly be- 
cause the higher-order diagrams in the current-current 
correlation function of Eq. often referred to as vertex 
corrections, are identically zero within this approach^ 

Thus, to evaluate the spin-diffusion coefficient we ap- 
ply the standard Matsubara formalism to Eq. using 
the above definition of the spin current: 



k 



dn{x) 



Ox 



(5) 



where k and v are the internal momentum and fre- 
quency of the "bubble" diagram, respectively. Here 
t>k = 2isin/c i a, A k (v) = — (l/7r)ImGk(^) is the electronic 
spectral function, n(v) = [eS u ~^/ T + is the Fermi 
function, and fj, is the chemical potential. 

The DMFT imposes a special form of the Green's func- 
tion in which the self-energy is k-independent and is de- 
fined from a self-consistency condition specified below. 
Within the DMFT parameters of the model are rescaled 
such that t = t^fz is finite as the dimensionality d — > oo, 
where z = 2d is the number of nearest neighbors^ In 
the following \plt is set to unity. Generally, the d = oo 
limit is well defined for Bethe and hypercubic lattice ge- 
ometries. While the semicircular electronic density of 
states (DOS) of the Bethe lattice is convenient for cal- 
culations, the Bethe lattice itself lacks the translational 
invariance and inversion symmetry implicitly used to ob- 
tain Eq. ©. Thus, it is important to determine whether 
the results on the Bethe lattice are equivalent to the re- 
sults on the hypercubic lattice, which is free from such 
deficiencies. Since many problems have been studied us- 
ing the Bethe lattice^ this comparison will have an even 
broader significance for the DMFT in general. 

We briefly sketch here the DMFT equations for the 
DE modeliiS* 1 ^ Since the self-energy is local, one can 
change G k (E) => G e (E) = [E — e — ^(E)}~\ and 
Sk ^ J^ e Po(e)j where e = £k, and Po(e) — \J2 — e 2 /n 



and Pq (e) = ey^(— e 2 ) / ^/n are the bare DOS's for Bethe 
and hypercubic lattices, respectively. The properties of 
the system are obtained from the local Green's function: 



de 



E-e-T,(E) 



(6) 



where the self-energy is defined from T,(E) = g^ l (E) — 
g^ 1 (E) 1 reminiscent of the Dyson equation, where 
g$ (E) is the "Weiss" function containing the dynamic 
influence of the environment on a given local site. The 
solution of the single-site problem provides a relation be- 
tween go(E) and S(_E). 1016 In the paramagnetic state 
and in the quasiclassical limit S > 1, such a relation is 
particularly simple^ £(i?) = ( J H S) 2 ga(E), which yields 



9(E) 



12(E) 



(JhS) 2 - Z(E) 2 



(7) 



Together with Eq. ©, this gives a self-consistent condi- 
tion for T,(E) or g(E). 

We now compare the Bethe and hypercubic solution 
for the single-particle properties. Fig. ^shows the evo- 
lution of the interacting DOS N(uS) — — (l/7r)Img(o;) for 
several values of JhS. The band splits as Jh increases 
and both lattice geometries exhibit the same qualitative 
behavior. Since the hypercubic DOS is expected to have 
in-gap states, one may ask whether the metal-insulator 
transition is well defined for a half-filled band. We find 
that the band splitting in the Bethe and hypercubic lat- 
tices happens at the same critical value (Jh&) c = 1/V2- 
At the transition, the imaginary part of the self-energy 
at lo = vanishes and the real part diverges. As a re- 
sult, N(w = 0) is exactly zero for Jh > J H and the 
metal-insulator transition for the half-filled band in the 
hypercubic geometry is well defined^ At small energies, 
N(u>) oc e~( JffS ) 1^ vanishes quite abruptly. 20 In con- 
trast, the "outer" tails of the upper and lower bands be- 
have very similar to the "bare" Gaussian form. Note that 
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FIG. 1: Interacting DOS N(ui) for several values of JhS for 
the Bethe (solid lines) and hypercubic (dashed lines) lattices. 



for JjjS > 1, the form of N(uj) for each subband is rather 
insensitive to the further increase of Jjj. 

We now rewrite Eq. © for the spin-diffusion coeffi- 
cient within the DMFT: 



0.6 



1 

2z 



— — — / dv 



dn(x) 



dx 



0_ 

db 



A{u) , (8) 



where b = -ImS(i^) and A(v) = J depQ(e)A £ (v) with 
Pq(e) being the "current" DOS defined from the trans- 
formation ^2y. sin 2 kia J depo(e)M- For the hypercu- 
bic lattice Po(s) — Pq(e)/2 while for the Bethe lattice 
pB(e) = (2 - e 2 )p B (e)/3. To obtain Eq. © we used the 
relation A e {v) 2 = [1/6— d/db] A E (v)/2n. Since the spin 
conductivity is proportional to the correlation function 
of two spin currents, each scaling as t ~ l/y/z-, the pref- 
actor in Eq. JSJ contains l/z. This means that dV 8 x is 
finite as d — > oo, similar to the particle conductivity^ 

In the manganites, the Curie temperature T c is much 
smaller than either the Hund's coupling or the band- 
width, in agreement with DMFT and Monte Carlo calcu- 
lations for the DE modeliifl* 2 ^ Therefore, all realistic tem- 
peratures are much smaller than the bandwidth T <C W 
and the derivative of the Fermi-function in the integrand 
of Eq. (jSJ) should be replaced by a 5-function at the 
chemical potential. Then, combining Eq. JSJl with the 
specific form of the DOS's for the Bethe and hypercubic 
lattices, one arrives at: 



V sX z \ a^(2-V(2-26 2 -/ 2 )-/.9' 



4 ,2-va-2& 2 )-2jy 



(9) 



where g' = Heg(v), g" = lmg(i>), and / = v — Re£(i/). 
Thus, the spin-diffusion coefficient is expressed through 
the local electronic Green's function and self-energy only. 
Fig. [2 presents T> s \ as a function of the electronic con- 
centration n for several JhS. As n varies from to 1, 
the chemical potential sweeps from u> = — oo to uj = in 
Fig. [2 The results are very similar in both geometries 
and become independent of JhS as JhS — > oo with the 
maximum located at n — 0.5. In the limit JhS 3> 1, 
Eq. lHU yields a numerical value for this maximum in 
V s xz/a 2 : 5/6tt(= 0.265) and 0.292 for the Bethe and 
hypercubic case, respectively. This demonstrates a close 
quantitative correspondence between the results in the 
Bethe and hypercubic lattices, which justifies the use of 
the former despite the concerns outlined earlier. An in- 
teresting feature appears in the results for the hypercubic 
lattice as n — > 0. Instead of vanishing, the spin-diffusion 
coefficient tends to a finite limit T> s xz/a 2 = 1/4tt(JhS) 2 , 
shown by circles in Fig. |2J 24 

In the high-temperature limit T > J H S > 1 Eq. © 
yields the result T> s x ~ l/T obtained previously using 
a Tchebycheff bounds (TB) formalism^ Numerically, 
(V° MFT /Vj B ) Bethe = 40\/2/9^ 3 / 2 = 1.13 agree very 
closely as well. 

It is interesting to analyze our results in the context 
of CMR materials. The superexchange (SE) interaction 
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FIG. 2: Vsxz/a 2 versus n for several JhS for the Bethe (solid 
lines) and hypercubic (dashed lines) lattices. 



is often discussed^ as necessary to correct the DE model. 
Within DMFT, the SE coupling J SE must scale as l/z. 
Combining this with the result of Ref. Q] one finds that 
V S S E oc l/z 3 / 2 is suppressed in comparison with the DE 
result T> s oc l/z. This also serves as a demonstration that 
the DE model cannot be simply reduced to an effective 
Heisenberg model. In a real d = 3 material the role of SE 
is further reduced by the smallness of jfj E ~ t 2 /JhS in 
comparison with the kinetic energy ~ xt. Therefore, the 
DE must dominate the spin diffusion and one can expect 
our results to be valid not only for the metallic part of 
the CMR phase diagram (x — 0.22. .0.5 for Lai^^Ca^Os, 
n = 1 — x in Fig. [2J but also for the ferromagnetic insu- 
lating phase (x < 0.22)^ Of course, in the limit x — > 
DE will diminish and the SE will dominate. Also, the 
critical scaling in the mixed phase < x < 0.12 (with a 
mixture of antiferromagnetic and ferromagnetic orders) 
can be strongly modified^ These cases require separate 
consideration. At x > 0.5 charge ordering prevents the 
DE mechanism from being operative. 

A systematic neutron scattering study of the metal- 
lic manganites recently focused on spin diffusion. 3 The 
width of the observed peak in <S(q, uj) centered at w = 
scales as Aq 2 , where A = 2T> S . Experimental re- 
sults taken at T ex P(x) ~ l.lT c (x), where T c (x) is the 
Curie temperature for a given hole concentration x, give 
A = 15 - 30 meVA 2 . To compare our results to ex- 
periments we use z = 6 and a — 3.87 A. But keep in 
mind that spin diffusion persists below the transition 
point and that, at least for the experimentally accessible 
wave- vectors, the correlation length saturates at about 20 
A. This implies that local magnetic correlations are sup- 
pressed by electronic inhomogeneities, which are prob- 
ably associated with the local charge ordering* 3 - So for 
simplicity, we take the susceptibility in the Curie form 
x ex P = S ,( S , + xyxr*, where S'(x) = S + (1 - x)/2 is 
an average on-site spin and T* — T exp (x) is known from 
Ref. 0- Fig. 01 shows the theoretical data for V s x/a 2 for 
JhS 3> 1 from Eq. @ as a function of x together with 
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FIG. 3: T> s x/a 2 as a function of x. Experimental values of 
T> s are from Fig. 5 of Ref. 0, x exp i s described in the text, 
theoretical results are from Eq. © for JhS 2> 1. 



experimental T> e s xp x exp / a 2 ■ We note here, that the theo- 
retical curves in Fig. are virtually independent of the 
actual value of JhS for JhS > 1. This figure demon- 
strates a remarkable agreement between the experimen- 
tal and theoretical results, which contains no fitting pa- 
rameters. Further improvement of the agreement can be 
sought, for example, from taking into account the second 
e g band which would effectively reduce J H /Wm& 



Since various other experiments also indicate the pres- 
ence of local inhomogeneities in CMR systems^ we pro- 
pose the following analysis. As q decreases, long-range 
magnetic correlations must dominate x(q) and the mag- 
netic correlation length must eventually exceed the size 
of the local polaronic distortions. So at a fixed tempera- 
ture close to T c , x e:cp (q) will increase and there will be a 
systematic decrease in the observed value oiT> e s xp ~ 1/x- 
This set of measurements would provide further infor- 
mation about magnetic correlations within the inhomo- 
geneities in CMR systems. 

In conclusion, we have presented a self-consistent study 
of the spin-diffusion in the double-exchange magnets 
within the framework of DMFT. This non-perturbative 
approach allows us to calculate the spin-diffusion coeffi- 
cient at any temperature down to a transition point. A 
good agreement with the experiments in the ferromag- 
netic CMR manganites and earlier work is found. Alto- 
gether, this provides a new insight into the dynamics of 
strongly-correlated itinerant magnets. 
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